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Abstract 

Pursuing a bottom-up approach to explore which flavor symmetry could serve as an 
explanation of the observed fermion masses and mixings, we discuss an extension 
of the standard model (SM) where the flavor structure for both quarks and leptons 
is determined by a spontaneously broken 84^ and the requirement that its particle 
content is embeddable simultaneously into the conventional 50(10) grand unifled 
theory (GUT) and a continuous flavor symmetry Gf like SO{3)f or SU{3)f. We 
explicitly provide the Yukawa and the Higgs sector of the model and show its viability 
in two numerical examples which arise as small deviations from rank one matrices. 
In the flrst case, the corresponding mass matrix is democratic and in the second one 
only its 2-3 block is non-vanishing. We demonstrate that the Higgs potential allows 
for the appropriate vacuum expectation value (VEV) configurations in both cases, if 
CP is conserved. For the first case, the chosen Yukawa couplings can be made natural 
by invoking an auxiliary Z2 symmetry. The numerical study we perform shows that 
the best-fit values for the lepton mixing angles 612 and O23 can be accommodated 
for normal neutrino mass hierarchy. The results for the quark mixing angles turn 
out to be too small. Furthermore the CP-violating phase 6 can only be reproduced 
correctly in one of the examples. The small mixing angle values are likely to be 
brought into the experimentally allowed ranges by including radiative corrections. 
Interestingly, due to the 5*4 symmetry the mass matrix of the right-handed neutrinos 
is proportional to the unit matrix. 
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1 Introduction 



The discovery of neutrino masses and attempts to understand the flavor puzzle have made 
it quite natural to expect the existence of an embedding product group for the SM such 
as SO{10) X Gf [1,2] at very high energies, where SO{10) acts as gauge and Gj as a flavor 
(family) symmetry. All fermions of one generation can then be unified into the spinor 
representation 16 of 5*0(10) and the three known generations are assigned to representa- 
tions of Gf. To specify Gf two important properties have to be fixed: Gf can be either 
abelian or non-abelian and it can be either continuous or discrete. Abelian symmetries are 
not able to explain the existence of more than one generation, since their representations 
are all one-dimensional. Therefore we consider such a choice of G/ to be less interesting 
than the others, although many models have obtained interesting results from U{1) fia- 
vor symmetries [3]. Most of the other models in which C/ is non-abelian [1] have some 
common features, like additional U{1) or Z„ factors, heavy vector-like fermions, elaborate 
mechanisms for the VEV alignment, to make them viable. Further unification might ex- 
plain these assumptions. We search for a simpler model by first constructing a low energy 
theory with the SM gauge group and a discrete non-abelian Gf and then showing possible 
embeddings of this theory into an SO{10) GUT and a continuous Gf. To really unify 
all three generations Gf should be either SO{3)f or SU{3)f. Our discrete symmetry will 
therefore be a subgroup of S0{3) or SU{3). We like to point out that none of the models 
mentioned in reference [1] with Gf being continuous follows this strategy. In cases where 
Gf is discrete [2], the symmetry is always broken at the GUT scale so that these models 
are also not comparable to our ansatz. Models in which the discrete non-abelian fiavor 
symmetry is only broken at low energies became very popular in the last few years [4-11]. 
However, these models can rarely be embedded into SO (10) x Gf. 

In working with discrete symmetries two issues are necessary to mention. First, the break- 
ing of a discrete global symmetry docs not lead to (unwanted) massless Goldstone bosons, 
unlike continuous symmetries. Second, if this breaking is only spontaneous, it might pro- 
duce domain walls which can be a serious problem [12]. It can be solved by either invoking 
low scale infiation or embedding the discrete symmetry into a continuous gauge symme- 
try [13]. A further issue which is only important, if the symmetry should be gauged, is the 
question of anomalies. Since we intend to embed the discrete symmetry into a continuous 
one at high energies, it is enough to make sure that the continuous one is anomaly-free. 
This can be checked by calculating the usual triangle diagrams. If the symmetry turns out 
to be anomalous, adding appropriate representations can solve this problem. 
Searching for an adequate discrete group we concentrate on the smallest subgroups of 
SO (3) or SU (3) which have at least one irreducible three-dimensional representation. is 
the smallest of these groups and has already been discussed extensively in the literature [5]. 
The second smallest group, sometimes called Ty, is rarely known, but its properties are 
in some sense similar to A4. Concerning the three-dimensional representations the same 
holds for the groups T' [10] (T' two- valued group of T and T isomorphic to A4) and (T/j 
isomorphic to x C2). In this paper, we focus on the group S4, which has a different 
group structure and is therefore worth exploring as a possible fiavor symmetry. This 
has already been done in several papers [6,7]. Pakvasa et al. [6] constructed an SM-hke 
model for quarks whose transformation properties under 5*4, however, do neither allow an 
embedding into SO{10) nor into S0{3)f {SU{3)f). An SO{10) model presented by one 
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of the authors (R.N.M.) and Lee [7] uses similar S'4 representations for the fermions and 
Higgs fields, but leads to the SMA solution to the solar neutrino problem which is ruled 
out ^. Recently, also Ma [6] discussed an S4 model, but he only considered the neutrino 
sector. Therefore there is still no working 5*4 model which addresses the flavor structure of 
all fermions in a unified manner and can be embedded into a GUT as well as a continuous 
fiavor symmetry. The aim of our study is to fill this gap. 

As noted, in this paper, we consider the non-abelian discrete fiavor symmetry S4 accom- 
panying the SM gauge group which naturally commutes with it. The gauge and the fiavor 
symmetry are assumed to be both broken only spontaneously at the electroweak scale. 
The particle content consists of the three known fermion generations together with three 
right-handed neutrinos and a set of six Higgs fields being SM doublets and transforming 
as + 2 + 3]^ under S4. The choice of the fermion content allows the embedding of our 
model into 5*0(10) x Gf. We observe that the number of parameters determining the mass 
matrices in our model, i.e. the Yukawa couplings and the VEVs of the Higgs fields, equals 
the number of observables, i.e. masses and mixing angles, in the CP-conserving case. If 
CP is violated, there are more parameters in the model than observables to fit. However, 
it is rather non-trivial to find parameter configurations which can reproduce all data, since 
5*4 constrains the mass structures as well as the Higgs potential. In our numerical study 
this turns out to be possible except for rather small deviations of the quark mixings which 
might be fixed by radiative corrections. 

The paper is organized as follows: Section 2 contains the group theory of 5*4, in section 
3 we present our model and argue why this is the minimal model with respect to the 
possible embeddings into 5*0(10) and Gf. Then, we display the structure of the mass 
matrices arising from the S4 invariant Yukawa couplings. After setting our conventions 
for the mixing matrices we give some numerical examples and comment on their viability. 
To complete the model, we calculate the Higgs potential and perform a restricted analysis 
of its possible minima. In section 6 we conclude and point to open questions in our 
model. Finally, the appendices contain Kronecker products, Clebsch Gordan coefficients 
and embeddings of S4 as well as some information on the Higgs potential and the selected 
minima. 

2 The S4 Group 

The group 5*4 is the permutation group of four distinct objects. It is isomorphic to the 
group O which is the symmetry group of a regular octahedron and so well-known in solid 
state physics. Its order is 24, i.e. it has 24 distinct elements. ^4 has five conjugate classes 
and therefore contains five irreducible representations which are all real. Among these are 
two one-dimensional ones, the identity (i.e. the representation being invariant under all 
transformations of S4, also called the symmetric representation) and the ant i- symmetric 
one (i.e. the one changing sign under odd permutations, also called alternating). In the 
following we will denote the identity one with 1^ and the anti-symmetric one with 3^2- 
There is one two-dimensional representation called 2 and two three-dimensional ones, 3^ 
and 32 ■ Out of these five irreducible representations only the two three-dimensional ones 

^The renormalization group extrapolations of neutrino masses in this model have not been studied and 
they might make a difference to the above conclusion since neutrinos are quasi-degenerate in this case. 
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Table 1: Character table of the group S4. (explanation see text) 



are faithful. Their characters Xi i-e- the traces of their representation matrices, are given in 
the character table, see Tabled There we use the following notations: Ci with ^ = 1, 5 
are the five classes of the group, °Cj is the order of the i^^ class, i.e. the number of distinct 
elements contained in this class, °hc^ is the order of the elements R in the class Ci, i.e. 
the smallest integer (> 0) for which the equation R = 1 holds. Furthermore the table 
contains one representative for each class given as product of the generators A and B 
of the group. From the generators A and B all other elements of 5*4 can be formed by 
multiplication. They ought to fulfill the following relations [14] : 

A^ = 1 , B^ = 1 and AB^A = B , ABA = BA^B . (1) 

We show one possible choice of generators in Appendix lA.ll Using them we calculate the 
Clebsch Gordan coefficients for all the Kronecker products. 

5*4 is the smallest group containing one-, two- and three-dimensional representations to- 
gether with the group T'. 

S4 can be embedded into 50(3) as well as in SU{3) (where it is isomorphic to the group 
A(24) [15]) and therefore gives the opportunity to embed our discrete flavor symmetry 
into a continuous one which is broken at a high energy scale. Possible embedding schemes 
are shown in Appendix IA.4I 

In our model the group 5*4 is broken completely at the electroweak scale, however this break- 
ing could also occur in two steps such that 5*4 breaks to one of its subgroups which is then 
completely broken. The non-abelian subgroups of 5*4 turn out to be already well-known 
as flavor symmetries: they are 5*3 (which is isomorphic to D^), and A4. Correlation 
tables containing the corresponding breaking sequences for the representations of 5*4 can 
be found in [16]. 
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3 The ^4 Model 



3.1 Particle Assignment 

In this subsection, we describe how to assign the fermions and Higgs bosons to different 5*4 
representations in such a way that the model can be embedded simuhaneously into 5*0 (10) 
and a continuous ffavor symmetry Gf. We argue that our choice which is displayed in Table 
El is unique in that sense. 

To embed the model into SO{10) all the fermion generations have to transform in the same 
way under 5*4 . Furthermore they have to transform either as trivial or as fundamental 
representation of the flavor group Gf, since all the other representations of these groups 
have a dimension larger than three. The only choice is then that they transform as 82 
under 5*4 apart from the trivial one where all the generations just form total singlets under 
the flavor group. An important point to note is that it is not possible to assign the three 
generations to 3^, since this representation cannot be identified with the fundamental 
representation of 50(3)/ or SU{3)f. Therefore our assignment is unique. 

In the next step, we have to choose the representations for the Higgs fields in order to 
give masses to the SM fermions. This choice depends on the desired Yukawa structure as 
well as the constraint to fill (a) certain representation(s) of the embedding group Gf. For 
fermions transforming as 82 under S4, the Higgs fields which can couple in an 5*4 invariant 
manner belong either to 1^, 2, 3^ or to 82, i.e. the only representation which cannot 
couple to form a total singlet under 5*4 is l2- Taking now for simplicity only the couplings 
which are symmetric in flavor space, i.e. lead to symmetric mass matrices, one is left with 
the representations 1^, 2 and 3^. Regarding the possible embeddings into 50(3)/ one 
recognizes that at least five Higgs fields transforming as 2 + 3]^ are needed and for an 
embedding into SU{3)f one needs six fields ~ ii + 2 + 3]^. Furthermore it turns out that 
the minimal version of five Higgs fields ~ 2 + 3]^ is not phenomenologically viable, since 
it leads to traceless symmetric mass matrices. Therefore the minimal setup of Higgs fields 
that we choose to get fermion masses contains six fields transforming as 1^+2 + 3^ under 
5*4. In the case of an embedding into 50(3)/ these representations are identified with 
1 + 5 and in the case of SU{3) / they are unified into the six-dimensional representation 
of 5f/(3)/ (see: Appendix rO|l . 

Neutrinos can also have Majorana masses apart from Dirac mass terms. Since the three 
right-handed neutrinos are also unified into one 32 under ^4 the only invariant mass term 
for the right-handed neutrinos is simply proportional to the unit matrix. The embedding 
of our model into 50(3)/ does not change this, however one has to keep in mind that in 
SU{3) f the irreducible three-dimensional representation with which 32 of 54 is identified 
is complex and therefore does not allow an invariant direct mass term for the right-handed 
neutrinos. Keeping our model as minimal as possible its embedding into 50(3)/ instead 
of SU (3) / is therefore preferred. Nevertheless, the inclusion of gauge singlets transforming 
as 6 under SU (3) / can give masses to the right-handed neutrinos, if the singlets acquire 
an appropriate VEV. Since the six- dimensional representation of SU{3) / contains a total 
singlet of 54 (see Appendix IA.4|) . this does not necessarily lead to the breaking of 54 at a 
high energy scale. Furthermore notice that the situation changes, if the model is embedded 
into a GUT like 50(10) at the same time. This will be discussed below. 
A non-trivial structure for the right-handed Majorana mass term requires the introduction 
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of gauge singlets which transform non-trivially under S4. To implement the canonical type 
I seesaw [17] the VEVs of such fields ought to be of the order of 10^^ GeV. In this case, the 
flavor symmetry S4 is broken at this high energy scale rather than the electroweak scale. 
Hence we discard this possibility. To keep the model as minimal as possible we also do not 
include SU{2)i Higgs triplets which could give rise to a type II seesaw mass for the hght 
neutrinos [18]. In some classes of 5*0(10) models this has become quite popular in order 
to get a relation between the light neutrino mass matrix and the difference of the down 
and charged lepton mass matrix [19]. Nevertheless fits using type I seesaw in SO {10) can 
also be found [20] . 

To embed the model as a whole into SO{10), the Higgs fields also have to be identified 
with certain 5*0(10) representations. In order to get a tree-level coupling they have to be 
part of either 10, 120 or 126 under 50(10). Since we fixed the Yukawa couplings to be 
symmetric in fiavor space the representation 120 drops out. Furthermore we observe that 
we need one 126 which transforms trivially under 5*4 for the right-handed Majorana mass 
term. We ought to choose its VEVs such that it does not give rise to a mass term for 
the left-handed neutrinos, since SU{2)l Higgs triplets arc absent in our low energy model. 
The minimal choice of fields would be: six Higgs fields transforming as (10, 1^ + 2 + 3]^) 
under (5*0(10), 5*4) and one 126 ~ 1-^ for the mass term of the right-handed neutrinos. The 
126 should also contribute to the Dirac mass term of the other fermions, since otherwise 
the masses of the down quarks and charged leptons were the same. In our opinion this 
is still not enough. Therefore at best one promotes each SM Higgs doublet to one 10 
and one 126, such that the 5*0(10) model has six ten- and six 126-dimensional Higgs 
representations. Among the 126s only the one which transforms trivially under 5*4 should 
develop a VEV at high energies and the rest only at the electroweak scale. 
To complete the model one needs at least one further Higgs representation, for example 
a 210. This representation together with the 5*4 invariant 126-dimensional representation 
should break 50(10) down to the SM with the Pati-Salam group as intermediate group. 
However it should not break the fiavor group 5*4 and hence has to be assigned to 1^ under 
5*4. With all these Higgs fields we believe that it is possible to make a viable high energy 
completion of our low energy model. 

In the numerical examples given below we fit the masses and mixing angles at the scale 
H which equals the W boson mass. In order to perform a fit at the GUT scale instead 
we would have to take into account the renormalization group running of all masses and 
couplings which will be complicated in a model with such a rich Higgs structure. 
Some issues still need to be discussed. Without constructing the 5*0(10) invariant Higgs 
potential it remains the question whether the advocated VEV configuration can be realized. 
Since our model contains several ten- and 126-dimensional representations and at least one 
2 10- dimensional one, there is a doublet- doublet splitting problem, i.e. one has to ensure 
that only six of the SM-like Higgs doublets have masses at the electroweak scale while 
the others acquire masses around the GUT scale. In the same manner one has to solve 
the well-known doublet-triplet splitting problem. In general separating the electroweak 
and the GUT scale will be difficult without supersymmetry (SUSY). Furthermore it has 
to be guaranteed that the gauge couplings unify at all being at the same time still in the 
perturbative regime at the GUT scale. 

Finally, the 5*0(10) x 5*4 model is embedded into 5*0(10) x Gf. As already mentioned, 
the three generations of fermions are identified with the fundamental representation of 
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Table 2: The particle content and its symmetry properties under S'4. 



Gf. The Higgs fields contained in the six ten-dimensional representations of 50(10) are 
unified into the two representations (10, 1) and (10, 5) for G/ = 50(3)/ and into (10, 6) 
for G/ = SU{3)f. The Higgs fields in the 126 s are treated in a similar way. Therefore the 
right-handed neutrinos acquire masses by coupling to ( 126 , 1) in case of Gf being 5*0(3)/ 
while they couple to (126,6), if G/ is SU{3)f. The 210 needed for the gauge symmetry 
breaking transforms trivially under O/ in the simplest case. Furthermore we need some 
50(10) gauge singlets breaking G/ in such a way that only ^4 remains. The smallest 
non-trivial possibility for G/ = 50(3)/ is 9 and for G/ = 5f/(3)/ 6, since they contain an 
54 singlet (see Appendix IA.4jl . But this we will not consider any further in this paper. 
To sum up, our model now contains three generations of fermions all transforming as 82 
and six SM-like Higgs fields transforming as + 2 + 3]^ under 54. In the next subsection, 
we show the arising Dirac and Majorana mass matrices. 

3.2 Fermion Masses 

The 54 invariant Yukawa couplings in our model are 

Cy = {Qi + Q2c' + Q3 1") ^0 + ai {Vs (Q2 - Qs t") ^1 + (-2 Qi + Q2 + Qs t") h) 
+ ((Q2 + Qs c^)li + {Qi + Q3 I2 + (Qi + Q2 u") Cs) 

+ a'^ (Qi + Q2S'' + Q3 b^) 4>o + af {Vs {Q2 - Q3 b") c^i + (-2 Qi + Q2S^ + Q3 b') 4>2) 

+ ai {{Q2 b'' + Q3 s=)6 + {Qi b' + Qs d^) 6 + (Qi + Q2 d') 6) 

+ (Ll + L2fi' + L3 t') + a\ (\/3 (L2 - ^3 r") + (-2 Li + L2 fi' + L3 t') ^2) 

+ ai ((L2 + L3 /i^)6 + (ii + L3 e') 6 + (Ll /i' + L2 e^) Cs) 

+ a"^ (Ll i/g^ + L2 z^;; + L3 i/^) ,^0 + (^ (L2 1^^^ - La z^^) 0i + (-2 Li z^,^ + L2 z^;^ + L3 z^^) ^2) 
+ ((L2 i^^: + L3 z.p6 + (Ll i.^ + L3 i.,^) 6 + (Ll z.^ + L2 z.,^) 6) (2) 

where the fields 0o,i,2 and ^1,2,3 are the conjugates of the fields 0o,i,2 and ^1,2,3 related by 
(f) = ecf)* with e being the 2-by-2 anti-symmetric matrix in SU{2)l space and the star 
denotes the complex conjugation. 
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They lead to the following mass matrices for i = u,d, e, u: 



^0 00 - 2 a\ 02 "2 ^3 "2 6 

M'={ aUs a^0o + ai(V3 0i + 02) «Ki | (3) 

"26 "Kl 0^0 00 + "1 (-v^0i + 02 

with the Higgs fields being replaced by their VEVs for the down quarks and the charged 
leptons and by the complex conjugate of their VEVs for the up quarks and the neutrinos. 
The sum of the VEVs has to be equal the electroweak scale, i.e. ^ |vEVjp ~ (174 GeV)^. 

i 

Note that the fields 0o,i,2 only appear in the diagonal entries. The contribution coming 
from the Higgs field 0o is proportional to the unit matrix, since 0o transforms trivially 
under S4. The fields 0i and 02 on the other hand are coupled in such a manner that 
their contribution is traceless. Finally, the fields which form a triplet under S4 only 
induce flavor-changing interactions, i.e. their contributions are encoded in the off-diagonal 
elements of the mass matrix Generally all the parameters in Eq.(j3)) can be complex. 

In the case of CP- conservation, we arrive at twelve real Yukawa couplings, five real VEVs 
for the Higgs doublet fields and the right-handed neutrino mass scale Mpt (see below). The 
sixth VEV is fixed by the electroweak scale. The 18 couplings and VEVs correspond to 
the twelve masses for the quarks, the charged leptons and the light neutrinos, the three 
CKM mixing angles and the three leptonic mixing angles whereof two have been measured. 
With CP- violation - either explicit or spontaneous - the number of parameters is increased 
such that it exceeds the number of observables. Nevertheless, it is not apparent that the 
mass structure restricted by the 5*4 symmetry allows one to fit all the data. Therefore, we 
perform a numerical study in the next section. 

It is interesting to note that assigning the fermion generations to 3]^ instead of 82 leads to 
exactly the same mass matrices, but does not allow an embedding into Gj without adding 
at least two further chiral generations to complete a representation of (see Appendix 

3. 



Next we discuss the neutrino sector. Since the right-handed neutrinos transform as 82 
under 5*4, their mass matrix Mnn is proportional to the unit matrix, i.e. Mjiji = Mnl. 
This means that the mass matrix for the light neutrinos arising from the type I seesaw has 
the form: 

= H^M^M''^ = H^iMr ■ (4) 

The last step is allowed, since all the Dirac mass matrices are symmetric by construction. 
The fact Mfui cx 1 indicates that the seesaw mechanism cannot be the sole origin of the 
difference between the quark and the lepton mixings in our model. As we will see below, 
in some cases there exists the possibility to impose an additional symmetry to maintain 
the diverse mixing patterns. 

With regard to leptogenesis, Mrj^ oc 1 is a viable starting point for the mechanism of 
resonant leptogenesis [21], since small radiative corrections can generate the small mass 
splittings which are needed. 

Finally , we want to indicate how the mass matrices change in case of a full 5*0 (10) model. 
First the number of Yukawa couplings is reduced from twelve {d)} with i = u^d^e^v and 
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j = 0, 1,2 to only six {ajji} {j = 0, 1,2 and R = 10 or R = 126 ), i.e. three couplings 
to (10, li), (10,2) and (10, S^) and three to (126, ll), (126, 2) and (126, 3i). At the 
same time the number of VEVs is in general increased. 

The additional embedding of S4 into a continuous flavor symmetry G / reduces the number 
of Yukawa couplings to its minimum. In case of G/ being SO{3)f the three Yukawa 
couplings {ajB,} for each R = 10 oi R = 126 can be expressed as two independent ones: 
aoR = C(R and aiR = 02 r = since the two- and three-dimensional representation under 
5*4 will be unified into one five- dimensional one. For Gf = SU (3) / we are left with only one 
Yukawa coupling for each i? = 10 or i? = 126, since the six-dimensional representation of 
SU (3) / breaks up into Ij +2-1-32 under 5*4. 

At the end, taking this setup does not reduce the number of parameters, but they can be 
eventually correlated such that predictions can be made. 



4 Phenomenology of the Mass Structures 

In this section we show that our model allows viable solutions. For this purpose, we 
introduce our conventions for the mixing matrices and give the experimental data. 



4.1 Conventions for the Mixing Matrices 

The Dirac mass matrices arise from the couphng: 
Vij Lf e (j) L^j 

for the downtype quarks (L = Q, L'^ = df^) and charged leptons (L = L, = e'^) and for 
the uptype ones (L = Q, L'^ = u^) and the neutrinos (L — L^IJ^ — u'^) : 

The mass matrices for the quarks are diagonalized by: 

UlM^'M'^Uu = diag(m^, ml, ml) , UlW^W^^Ua = diag(m^, m^ ml) 

where C/„ and are the unitary matrices transforming the left-handed up and down 
quarks to their mass eigenstates. The CKM mixing matrix is given by: 

VcKM ^U^U^ . 

The standard parameterization for Vckm is [22]: 

/ C12C13 S12C13 siae"*'^ 



Vr 



CKM 



-S12C23 - C12S23S13 e'" C12C23 - S12S23S13 e'" S23C13 
V S12S23 - C12C23S13 e*'' -C12S23 - S12C23S13 e^^ C23C13 



where we use the abbreviations = sin(^jj) and = cos(^jj). The angles are restricted 

to lie in the first quadrant and 5 can take any value between and 2 vr. 

Similarly, in the leptonic sector the mass matrix for the charged leptons fulfills the relation: 

UlM^M^^'Ui = diag(m^, mj,m^) 
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with the unitary matrix Ui transforming the left-handed charged leptons to their mass 
eigenstates. The hght neutrino mass matrix coming from the type I seesaw is a complex 
symmetric matrix. Hence it can be diagonalized by Uiy-. 

UlM.U* = diag(mi, m2, mg) , UIM.MIU, = diag(m2, ml, ml) 

with the mass eigenvalues being positive definite. The Umns niatrix is defined as 

3 

i^aL = ^ Umns ^iL for a ^ e, fx, r and i = 1, 2, 3 

i=l 

where the denote the flavor and the the mass eigenstates. 
Therefore Umns is expressed as 

Umns = Uf U* . 

In the case of Majorana neutrinos the Umns matrix can be factorized in a unitary matrix 
which is parameterized in the same way as the CKM matrix Vckm and a diagonal matrix 
containing the two Majorana phases (pi and ip2- 

Umns = Vckm ■diagie''^\e''^\l) ■ 

Both Majorana phases are taken to fulfill <(/?!, </72 < tt. 

4.2 Experimental Data 

The quark masses at the scale fi ^ Mw are given by [23]: 

mu = 2.2 MeV , mc = 0.81 GeV , m* = 170 GeV , 

ma = 4.4 MeV , = 80 MeV , mb = 3.1 GeV . 
The mixing angles and the CP-phase measured in tree- level processes only are [22]: 

si2 = 0.2243±0.0016 , S23 = 0.0413±0.0015 , S13 = 0.0037±0.0005 , 5 = 1.05±0.24 . 

They are almost independent of the scale /i at low energies. To quantify the CP-violation 
one can introduce the Jarlskog invariant J'cp [24]: 

= (2.88 ± 0.33) X 10-^ . 

The charged lepton masses at // Mw are: 

me = 511 keV , m^ = 106 MeV , = 1.78 GeV . 

In the neutrino sector only the two mass squared differences measured in atmospheric and 
solar neutrino experiments are known [25]: 

Amj^ = ml-m\ = (7.91°;^) x IQ-^eV^^ , |Am^^| = \ml-m\\ = (2.21°;^) x lO'^eV^ . 

The leptonic mixing angles are constrained by experiments: 

< 0.031 , — 0-3iao5 ' ^23 ~ 0-^iai2 ■ 

All values observed in neutrino oscillations are given at 2(7 level. The three possible 
CP-phases 5, </7i and </72 in the leptonic sector have not been measured till today. 
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4.3 Numerical Examples 



In this section, we present two numerical examples which can reconcile the data apart from 
the fact that some of the quark mixings turn out to be smaller than the central values. 
They correspond to perturbations around two different rank one matrices for the quarks 
and charged leptons: 



i.e. the matrix Aii is the democratic mass matrix, which has been discussed several times 
in the literature [26], and only has a non- vanishing 2-3 block. Matrices of the form 
of are often used as approximation for the light neutrino mass matrix to produce 
maximal atmospheric mixing [27]. Both matrices are able to explain the strong hierarchies 
observed in the quark and the charged lepton sector. For the neutrinos it is necessary 
that their mass matrix significantly deviates from the forms of A^i and to satisfy 
the restrictions on the lepton mixing matrix. In particular in our two numerical examples 
the neutrinos dominantly couple to the Higgs fields 0i,2, i-e. their mass matrix has large 
entries on its diagonal. 

In the following we show that certain VEV configurations which reflect our flavor symmetry 
together with some fine-tuning of the Yukawa couplings allow the matrices Aii and 
to result from the general form given in Eq.Q. 

First, we consider the matrix Aii. In the CP-conserving case, the Higgs potential has a 
minimum at which the VEVs of the fields are equal, the VEVs of 0i,2 are zero and (po has 
a non-vanishing VEV. The VEVs of and the field 0o are not necessarily equal, but we 
make the natural assumption that they are almost the same. With this and the constraint 
that the Yukawa couplings and a\ should be the same, the Dirac mass matrices 
have the democratic form. We then perturb around this known minimum of the Higgs 
potential and also allow changes in the Yukawa couplings a* in order to get mass matrices 
whose masses and mixing parameters agree with the observed ones. For concreteness: 



«0 


= 0.651341 -0.00001 i , < = -0.0058575 + 0.001286 i , = 0.651341 , 


< 


= 0.011598 + 0.000219 1, 


af = -0.00299585 - 0.00098905 i , 


«2 


= 0.011708 -0.0002841, 




«0 


= 0.0071585 - 2.0- 10~M 


, al = -0.000375685 + 0.00331405 1,^2 = 0.0065433 


«0 


= 0.15224 + 0.09061, < 


= 1.04426 - 0.018245 1,^2 = 0.08364 + 0.04734 i , 


(0o) 


= (87 + 0.01776 1) GeV, 


= (-2.5912 - 14.1982 i) GeV , 


(02) 


= (1.27- 7.6236 1) GeV, 


(6) = (86.8607 + 1.14585 i) GeV , 


(6) 


= (87- 0.11162 1) GeV, 


(^3) = (87 + 0.88384 i) GeV , 


Mr 


= 4.3 X 10^^ GeV . 





The hierarchies among the Yukawa couplings ~ ^ <y\ for i = u,d,e and ^ 
^2 ~ c^o can nicely be explained by an approximate auxiliary Z2 under which only the 
right-handed neutrinos and the Higgs fields (pi and 02 transform: 



-i^t and 01,2 -01,2 (6) 




(5) 
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whereas the rest remains invariant. The structure of the resulting mass matrices is then 
democratic for the quarks and the charged leptons while it is dominated by the (2, 2) entry 
for the light neutrinos. The numerical values given here lead to: 

niu = 2.2 MeV, m, = 0.814 GcV, = 169.94 GcV , 

nid = 4.46 MeV , m, = 81.2 MeV , rub = 3.05 GcV , 

me = 514keV, m^ = 105.5MeV, m^ = 1.76GcV, 

Aml-^ = 2.3 X 10"^ eV^ , Aml^ = 7.89 x 10"^ eV^ . 

The sum of the three light neutrino masses is rrii = 0.0614 eV, i.e. the mass spectrum 
is strongly hierarchical with the smallest mass mi 0.0038 eV. This is well below the mass 
bounds known from cosmology: Ylii'^i < (0.42 . . . 1.8) eV. This upper bound depends on 
whether the measurements of the Lyman a spectrum are included or not [28] . 
The quark mixing angles turn out to be 

si2 = 0.2238 , si3 = 0.003694 , S23 = 0.02831 . 

Unfortunately. §23 is too small. This may be compensated by radiative corrections. The 
CP-phase 5 is 1.174 radian and therefore near the upper bound 1.29 radian. The Jarlskog 
invariant has the value Jcp = 2.1 x 10~^. More interesting to see is that the neutrino 
mixing angles can be accommodated: 

s?2 = 0.3 , = 0.49946 , 

and we predict \U'^^g\ to be 0.06616. This is within the reach of the next generation 
experiments [29]. Furthermore, the three leptonic CP-phases are: 

S = 0.9983 , ipi = 1.579 , ip2 = 1.336 in radian. 

Calculating for completeness the magnitudes of \mee\ and m^ which can be extracted from 
neutrinoless double beta decay and beta decay, respectively, one finds: 

3 / 3 \ 1/2 

\mee\ = I Yl {UmnsT = 0.0054eV and mp = Ij^ I^mat^I'^? = 0.0069eV 

i=l \i=l / 

These values are at least two orders of magnitude below the current experimental bounds 
which are |mee| < 0.9 eV [30] and m^ < 2.2 cV [31]. They are also below the limits of the 
experiments planned for the next years [32,33]. The main reason for this is the strong 
hierarchy in the neutrino mass spectrum. 

If we take the second mass matrix A^2 as starting point, we find the following numerical 
example: 

= 0.56672 - 0.00001 i , < = 0.2833 + 0.00001 i , = 0.85175 + 0.00608 i , 
= 0.01028 - 3 • 10^^ i , < = 0.005145 - 0.000021 i , = 0.015597 + 0.000789 i , 
a^o = 0.0059333 - 9 • 10~^ i , = 0.0029676 - 10"^ i , = 0.0088571 + 0.0010664 i , 
= 0.00333 + 0.0281, < = 0.51567 - 0.01616 i , = -0.01572 + 0.89947 i , 
(00) = (99.9974 + 0.0026 i) GeV , (0i) = (-2.67385 - 6.79332 i) GeV , 
(02) = (99.9923 - 0.02867 i) GeV , (^i) = (99.9907 - 0.266 i) GeV , 
(^2) = (0.0058 - 0.38774 i) GeV , (^3) = (0.04332 - 0.14671 i) GeV , 
Mr = 6.3 X 10^^ GeV . 
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Here, one can clearly see that the Yukawa couplings have to fulfill the relation: : a\ : 
ttg ~ 2 : 1 : 3 for i = u,d,e to produce a matrix with a dominant 2-3 block. The resulting 
masses and mass squared differences are given by: 

m„ = 2.4 McV , = 0.812 GeV , rrit = 170.24 GeV , 

rud = 4.4 McV, m^^SOMeV, m6 = 3.10GcV, 

me = 512.6 keV , = 106 MeV , = 1.78 GeV , 

Aml^ = 2.4 X 10"^ eV^ , Aml^ - 7.59 x 10"^ eV^ . 

The mass spectrum for the three hght neutrinos is normally ordered and degenerate with 
mi ~ 0.1682eV and X^i"^* — 0.512 eV which coincides with the upper bound of the 
cosmological measurements, if the Lyman a data is also taken into account. For the 
quarks all three mixing angles are about 10% too small: 

si2 = 0.2128 , si3 = 0.0038 , S23 = 0.0389 . 

Again, this has to be compensated by radiative corrections. The CP-phase 6 turns out 
to be the more severe problem, since S fa 0.386 radian, i.e. the CP-violation generated 
by this setup is more than a factor of two too small and so is the Jarlskog invariant 
J'cp = 1.15 X 10~^. The two measured mixing angles in the lepton sector can be adjusted 
to their currently given best-fit values, i.e. = 0.306 and S23 = 0.506. The third mixing 
angle 6*13 and the three leptonic CP-phases turn out to be: 



Again, all phases are given in radian. \mee\ and mp have almost the same value 0.168 eV, 

since the neutrino masses are nearly degenerate and all phases are approximately tt. 

is near the limit which will be reached by the KATRIN experiment [32]. Due to the 

degeneracy of the neutrino mass spectrum also |mee| [33] and the sum of the neutrino 
masses [34] can be measured by the up-coming experiments in the next five to ten years. 

The matrix structure is very similar to the one of for the quarks as well as for the 
charged leptons while the matrix for the light neutrinos has approximately the form: 



In this case one also finds that the invoked VEV configuration = (^2) = (Cs) = 
together with the other VEVs being nearly the same is an allowed minimum of the Higgs 
potential, if CP is conserved. 

In the two examples shown here we have taken all parameters to be complex, i.e. we 
have assumed explicit CP-violation and have not made use of possible field re-definitions 
to absorb some of the complex phases. We have done so in order to keep as many (free) 
parameters as possible in the numerical fit procedure. For example, this does not mean 
that there is no viable solution in the case of spontaneous CP-violation. 



5^3 = 0.0034 {\U\ 



re3 I _ 
MNS\ — 



0.0584) , S = 3.032 , cpi = 3.102 , cp2 = 3.081 , 




(7) 
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5 The Higgs Potential and its Possible Minima 



Finally, we discuss the S4 invariant Higgs potential V and show that there exist two 
CP-conserving minima from which the VEV configurations assumed above can arise as 
"perturbations" . 



V 



(8) 



+ Ao + Ai {4h + 4'p2? + A2 i4>\'p2 - 4<Pif 



+ A3 



+ ai {4o'f>o){4l'Pl + 4'^2) + {(T2 [(4</'l)^ + i4o'p2f] + h-C.} 



+ a2 



[|4</.l|2 + 14,^21'] + {^3 [i4o^l)i4l<t>2 + 4'^l) + i4o^2){4lh - 42<f>2)] + h.C.} 



+ ^1 (Y. j + ^2 [3 (ek2 - d^sf + (-2 ^ki + + d^3? 

+ 4 [(^ks + cla? + (ck3 + cki)' + (d6 + cki)'^ 



+ Ai 



(eks - ^k2)' + (ek3 - eki)' + (^k2 - eki)' 



+ n(0S<^o)(E^l^O+^2 {Y.4<p, 



+ T3 [\/3(</>i</>2 + 4'Ai)(ek2 - cks) + ('/'I'^i - 4'/'2)(-2eki +ck2 +?k3) 
+ {m [4 (4a)' + (V3 46 + 46)' + (V3 46 - 46)'] + h-c} 



-l-h.c 



■} 



+ ki [4 |46P + |x/346 + 461' + |x/346 - 461' 

+ {k2 [4 (4^1)2 + (^/346 - 46)' + (^^46 + 46)' 

+ k2 [4 1461' + 1 V3 46 - 461' + 1 V3 46 + 461'" 
+ [2 (46)(ck3 + d6) - (V3 46 + 46)(ek3 + cki) + (^46 - 46)(ek2 + eki) 

+ h.c.} 

+ [2 (46)(?k2 - eks) + (^^46 - 46)(6k3 - eki) - (^46 + 46)(eki - 6^2)' 

+ h.c.} 



{^^5 [(46)' + (46)' + (46)' 



+ h.c. > + K5 



1461' + 1461' + 1461' 



+ 

+ {'^e [(46)(?k3 + ek2) + (46)(6k3 + eki) + (46)(6k2 + eki)] + ^.c.} 

+ [u, [v^(40i)(^t^2-^k3) + (4'^2)(-26ki+ck2+ek3)]+h.c.} 

+ {c^2 [2(46)(46) - (46)(V3 46 + 46) + (46)(v^46 - 46)] + h.c.} 

+ {a;3 [2(^I</>o)(46) - (e|</'o)(V346 + 46) + (el0o)(V346 - 46)] + h.c.} 
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where parameters with +h.c. in curly brackets are in general complex (for example the 
last parameter CJ3) and the rest is real. The Higgs potential has 30 parameters in total. 
27 of them are quartic couplings out of which 11 are complex. Interesting to notice, the 
potential is invariant under the following transformation: 



01 -01 and ^2 ^ ^3 

and the fields 0o, 02 and ^1 remain unchanged. 
We can parameterize all possible real VEVs as: 



(9) 



(10) 



(0o) = ^"0 , (0i) = u cos{a) , (02) = u sm{a) , 

(^1) = V cos(/3) , (^2) = V sin(/3) cos(7) , (^3) = v sin(/3) sin(7) . 

The potential at the minimum has then the following form: 

Vmin = -lA - ulu^ - l4, + Ao + (Ai + A3) (11) 

+ {\{ + 4 A| sm2(/3)) v^ + \l (2-3 sm^{P)f + 3 sin^(/3) 008^(27) 

- A| (3 + cos^{2 7)) sin'^(/3) + {ai + 2 Re(cr2) + ^2) vlu^ + 2 Rc{a^) sin(3a) vq 
+ (n + 2Re(«;5) + K5) vlv'^ + (4Re(Ki + K2) + 2 (ki + R2) + T2) v? 
+ (2 Re(«;i - K2) + ki - K2 + ^3) [- cos(2 a) (2-3 sin^(/3)) 

+ \/3 sin(2a) sin2(/3) cos(2 7) v'^ 

+ 3 Re(K6) sin(/3) sin(2 /?) sm(2 7) vq 



+ 2 Re(a;2 + cos — toi) sin(a ) (2 - 3 sin2(/3)) - V3 cos(a) sm'^{P) 003(27) 



UVqV 



Note that the couplings A2, A4, , do not appear in Vmin- One can deduce the following 
VEV conditions: 
dV ■ 



da 



2u 



V^Vo 



+ 



dVrr, 



d/3 



cos(a) (2-3 sin2(/5)) + Vs sin(«) sm^{j3) cos(2 7)j yi (12a) 

uv^ (sm{2a) (2-3 sin^ip)) + V3cos{2a) sin^ip) cos(2 7)) y2 
+ 3vou^ cos(3 a) Re((T3)] 

= v'^ —2V3uvo ^-\/3 sin(Q;) + cos(q;) cos(2 7)) sin(2/3)yi 

cos(2 a) + sin(2 a) cos(2 7)) sin(2/3)7/2 

+ 2 (sin(4 p) + sin^ip) sin(2 /3) sin2(2 7)) 7/3 
3 

+ -Wo {S sin(3 /9) — sin(/3)) sin(2 7) Re(«;6) 



(12b) 



97 



= 2^^ sin^(/?) 2\Piuvo cos(q;) sin(2 7)^/1 

— -s/Sii^ sin(2Q;) sin(2 7)y2 
+ sin^(/3) sin(4 7)^/3 

+ ^vvq cos{P) cos(2 7) Re(«;6)] 



(12c) 
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with Ui being defined as: 




(13a) 
(13b) 
(13c) 



All three equations Eq. ()12|l have to be equal zero. A more restrictive requirement is that all 
the terms should vanish separately. The Re(K6) and Re(cr3) are parameters of the Higgs 
potential and should not be constrained to vanish in order to avoid accidental symmetries 
arising in the potential. Therefore their coefficients should vanish separately. This poses 
restrictions on the angles a, j3 and 7 as well as on the moduli f , u and vq. Obviously, one 
of the solutions is given by u = (a is then no longer a variable), (3 = arccos(l/\/3) and 
7 = 7r/4 , i.e. the VEVs of the fields 0i,2 vanish, the VEVs of the fields are equal and 
00 has an in general non- vanishing VEV vq. Assuming that all quartic couplings are of the 
same order and all mass parameters have the same order, it is natural that the VEV for 
the fields is of the order of Vq. This means only a slight parameter tuning is necessary 
to achieve the equivalence of these VEVs as is needed for the zeroth order approximation 
of the fermion masses in our first numerical example. It is noteworthy that the number of 
(massless) Goldstone bosons is increased by two, if one additionally sets the parameters 
Hi, Re((T3) and Re(K6) to zero, see Appendix IB. 21 I.e. the restrictive requirement that 
none of them vanishes turns out to be sufficient to avoid further Goldstone bosons. In the 
numerical study it is pointed out that an auxiliary Z2 symmetry can explain the required 
Yukawa couplings. This Z2, if also valid in the Higgs sector (and therefore in the whole 
Lagrangian), restricts the quartic couplings. It enforces as, /ta, K4, cji, uj2 and uj^ to vanish. 
As far as we can see this does not create an accidental symmetry in the Higgs potential 
and so does not change the discussion. 

A similar analysis can be done for our second numerical example which enforces the VEVs 
of 00; 02 and ^1 to be equal and the other VEVs to vanish in the zeroth approximation. 
This corresponds to a = | and /? = (7 is then irrelevant). One sees that also in this 
case the coefficients of the parameters Re(K6) and Re((j3) vanish such that the VEV 
conditions Eq.lfT^ can be fulfilled. Note that the values of Wo; u and v are not constrained 
to be equal, but again it is plausible that they are nearly the same, if the mass parameters 
IXi as well as the quartic couplings of the Higgs potential are chosen to be of similar size. 
Interestingly, setting the parameters Ui, Re(/€6) and Re(cr3) to zero increases the symmetry 
of the potential also at this minimum and leads to three further Goldstone bosons (see 
Appendix IB.3|) which is actually one more than in the case above. By inspecting the Higgs 
potential V one finds that at least for real VEVs and parameters of the potential requiring 
Hi = 0, Re(cr3) = and Re{nQ) = causes an accidental S0{2)acc symmetry under which 
(01, 02) forms a doublet and the other fields remain invariant and an accidental S0{3)acc 
under which (^1,^2,^3)^ transforms as triplet and the fields 0o, 0i, 02 trivially. The 
VEV configuration with vanishing (0i,2) then only breaks S0{3)acc and not S0{2)acc and 
therefore gives rise to two Goldstone bosons whereas the configuration (0o) 7^ 0, (02) 7^ 
and (^1) 7^ breaks both accidental symmetries resulting in three Goldstone bosons. 
In the limiting case that all mass parameters fii of the potential are equal, these two 
minima are exactly degenerate (along with many others), since the value of the potential 
at the minima is — | (/i^ Vq + fil t>^) and — | (/i^ Vq + filu^ + /ig v"^), respectively. 
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Further investigation of the minima with CP-violation which lead to realistic masses and 
mixing parameters is beyond the scope of the paper, but it is plausible that these minima 
can be formed through small deformations starting with CP-conserving minima, as done 
here. 

We did not perform any checks of the stability of each minimum and the potential as a 
whole, since the number of parameters (~ 30) makes us confident that there exists at least 
one point in the parameter space for each minimum where it fulfills together with the 
potential all the stability criteria. Furthermore we did not address the question whether 
the minimum is a local or global one, since this might also only depend on an appropriate 
choice of the parameters. 

6 Conclusion and Outlook 

To conclude, we have presented a low energy model based on the SM gauge group aug- 
mented with the flavor symmetry S4. In contrast to other flavor models we used the 
requirement to embed our model into a GUT like 5*0(10) and at the same time into a 
continuous flavor group hke SO{3)f or SU{3)f as guideline for the choice of the transfor- 
mation properties of fermion and Higgs fields under S4. The resulting model is minimal 
in that sense. 

Since the structure of the mass matrices is determined by 6*4, it is not obvious whether we 
can accommodate all observed masses and mixing angles, even though the model contains 
as many parameters as observables needed to be fixed in the CP-conserving case. 
To check this, we explore two cases which are perturbations around two different rank one 
mass textures for quarks and charged leptons that can be maintained for two choices of 
ground states of the theory together with some tuning of the Yukawa couplings. The first 
is the so-called democratic mass matrix and the second one only has a non-vanishing 2-3 
block. We give numerical examples for each that are able to fit the known fermion masses 
and mixing angles in the quark and lepton sector up to rather small deviations. We believe 
that invoking radiative corrections will lead to full accordance with the experimental data. 
The difference between the mixings of quarks and leptons crucially depends on the fact that 
the form of the mass matrix of the light neutrinos differs strongly from the one of the quarks 
and charged leptons. In our first example an auxiliary Z2 can help to explain this difference 
and in the second one the parameters have to be fine-tuned. Taking the auxiliary Z2 as an 
exact symmetry of the theory prevents the model from being embedded into SO (10), since 
the right-handed neutrinos transform differently from the other fermions. However, one 
can still promote our model to an 5*^7(5) GUT. As the Higgs fields 0i^2 transform under 
Z2 whereas ^1,2,3 remain invariant, the Z2 is not compatible with any embedding of our S4 
fiavor symmetry into SO{3)f and SU{3)f without adding further fields. 
The right-handed neutrinos are degenerate at tree-level and even more their mass matrix 
is proportional to the unit matrix. Therefore the large leptonic mixing angles have to be 
encoded in the structure of the Dirac mass matrices for the neutrinos and charged leptons. 
The VEV configurations we used in our numerical examples can only be analyzed in the CP- 
conserving hmit, since the Higgs potential turns out to be quite complicated. Nevertheless 
these are determined by our flavor symmetry. In contrast to this, the values of the Yukawa 
couplings are not fixed by 5*4. The question why the top quark is 36 times heavier than 
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the bottom quark remains unanswered, but can possibly be explained, if our model is 
promoted to 5*0(10) x Gj. 

Throughout this paper we have not been concerned with the question how to guarantee 
that the Higgs spectrum just contains one light uncharged (scalar) Higgs inducing only 
flavor diagonal interactions like the Higgs in the SM while the rest is heavier. Connected 
to this is the problem of flavor changing neutral currents and lepton flavor violations which 
usually arise in models with more than one SM-like Higgs doublet. Typically these effects 
are negligible, if the masses of the flavor changing Higgs flelds are above a few TeV. 
Systematic calculations are difficult, since the Higgs mass spectrum cannot be evaluated 
in general cases. 

Finally, we want to comment on the possibility of supersymmetrizing our model. Intro- 
ducing supersymmetry apart from its salient feature to solve the hierarchy problem if it is 
broken at low energies technically leads to a severe simplification of the Higgs potential, 
since then all the quartic terms are determined by the D-terms. The danger lies in the 
fact that this generally leads to large accidental global symmetries in the potential which 
consequently lead to a number of unwanted Goldstone bosons. Two ways of treating this 
problem can be found in the literature: a.) breaking the discrete and hence also the acci- 
dental symmetries by the soft SUSY breaking terms (for example: [35]) or b.) introduce 
gauge singlets whose couplings are invariant under the discrete symmetry, but break the 
accidental ones (for example: [36]). Obviously, the whole situation can change in a grand 
unified model, since then the Higgs doublet fields can belong to various representations 
of the GUT which have different invariant couplings (see for example [37] for a SUSY 
50(10) model). Since supersymmetric potentials are restricted to be positive by construc- 
tion, checks of their stability are easier than for non-supersymmetric ones. These issues 
are currently under study. 
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In this appendix we display the representation matrices, Kronecker products and Clebsch Gordan 
coefficients to calculate all the terms being invariant under the group 5*4. 

A.l Representation Matrices 

The representation matrices fulfilHng Eq.(^ can be chosen as: 



A Appendix: Group Theory of S4 




and 




for 2 



A = 




and B = 
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and 

/ 1 

and B = 1 I for 82 . 

\ 1 

These matrices can be found in [14]. 

A. 2 Kronecker Products 

The Kronecker products can be calculated from the above given character table [38]. 

liXlj = l(i+j)„,od2+i Viandj 
2x1, = 2 Vi ^i^^ = ^1+^2 Vi 

0^1 o vy ^ 3ix32 = 12 + 2 + 31 + 32 

3i X Ij = 3(i+j) 2 +1 V 1 and J 

[2x2] = li + 2, {2x2} =12 and [Sj x 3i] = li + 2 + 3i , {3ix3i}=32 Vi 

where we introduced the notation [fj, x fj] for the symmetric and {fj, x /x} for the anti-symmetric 

part of the product ^ x 11. 

Note that 1/ x = x 1/ for all representations jj, and v. 

A. 3 Clebsch Gordan Coefficients 

The Clebsch Gordan coefficients can be calculated [39] with the given representation matrices for 





b[ 




. b'2 


bs J 


\ b's 









5)1 


c'2 


1 ~32 









Since we choose all the representation matrices to be real, it also holds: 

~ ll , ~ I2 , ) ~ 2 , 1^ ) ~ 3i and ( cl | ~ 82 

The Clebsch Gordan coefficients for the one-dimensional representations are trivial: 

as well as the products 1^ x /x of any representation /it with the total singlet Ij : 

f Aai 
V Aa2 





^1 , -^C2 ~02 
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And here are the ones for I2 x fJ- oi any representation /x: 



-Ba2 
Bai 



Bh\ / Bci 

-32 , Bc2 

Bh \Bc3 



The Clebsch Gordan coefficients for fi x ji have the form: 



for 2 



aia'i + 02^2 ~ il 



-aia'2 + a2fl'i ~ I2 



for 3i 

3=1 



aia'i - a2a'2 y " - f 62^3 + ^3^2 \ / ^3^2 

61^3 + ^3^1 ~ 3i , hbi, ~ b3b[ I ~ 32 




for 32 



E ~ il 



C2C3 + C3C2 

C1C3 + C3C1 I ~ 3i , I - csc'i I ~ 32 




C1C2 + C2C']^ 

Note here that the parts belonging to the symmetric part of the product n x fi are symmetric 
under the interchange of unprimed and primed whereas the ones belonging to the anti-symmetric 
part change sign, i.e. are anti-symmetric. 

Note also that for our choice of generators the Clebsch Gordan coefficients for 3^ x and 
^2 X 32 turn out to be the same. For the coupling 2 x 3^ the Clebsch Gordan coefficients are 

0261 \ / aici 

-^(\/3ai62 + 02^2) )~3i 1 ^{V3a2C2 - aiC2) ) ~ 3^ 

l(V3a,bs-a2bs) J and for 2 x 82 V -^(V3a2C3 + 0^03) 

aibi \ / 02C1 

i(\/3a262 - 0162) ~32 -i(\/3aiC2 + a2C2) |~32. 

-i(A/3a263 + ai&s) / \ ^iV^aiCa - a2Cs) 

And for 3^ x 82 one finds the following combinations: 

3 

E bjCj - I2 



^(26ici - 62C2 - 63C3) 
^{b2C2 - 63C3) 
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63C2 - 62C3 

6lC3 - 63C1 
62C1 - 61 C2 



&2C3 + 63C2 
^'1C3 + 63C1 
blC2 + ?'2Cl 



32 



A.4 Embeddings of ^4 into SO(S) and SU(S) 

We only display the resolution of the smallest representations of 50(3) {SU{3)) into irreducible 
ones of S4. 



S0{3) 
1 
3 
5 
7 
9 



S4 
ll 
32 

2 + 3i 

12 + 31+32 
ll + 2 + 3i + 32 



SU{3) 
1 
3 
6 
8 
10 



5*4 

ll 
32 

11 +2 + 3i 
2 + 3i + 32 

12 + 3i +2 32 



The first table can be found in [14] and the second one can be calculated with the formulae given 
in [40]. 



B Appendix: Minimization of the Higgs Potential 
B.l Remaining VEV Conditions 

The derivatives and have the following form: 



dVrr, 



dv 



dVrr, 



dvn 



dVrr, 



du 



= 2v (-/xi + 2(Af + 4A|)u2) +2vug(2Re(K5) + K5 + Ti) 

+ 2vu^(4Re(«;i + «;2) + 2(Ki + K2) + T2) + 9 vqv Re(/C6) sin(/3) sin(2/3) sin(27) 



+ Auvqv 



yi 



+ 2u^v 



y2 



sm{a) (2-3 sin^ip)) - Vs cos{a) sin^iP) 003(27) 
- cos (2a) (2-3 sin2(/3)) + V3 sin(2 a) sin (P) cos(27) 

+ 4t.3 [4 sm\(3) - (3 + cos2(2 7)) sin^(/?)] 7/3 

= 2 (-/x? + 2 Ao Uo) + 2 vo v? [ox + 2 Re(cr2) + 0^2) + 2 Re(c73) sin(3 a) 
+ 2i;o'u^(2Re(/«5) + K5 + Ti) + 3i;^Re(AC6) sin(^) sin(2/3) sin(27) 



+ luv^ 



sin (a) (2-3 sin2(/3)) - cos (a) sin2(/3) 003(27) 



2u (-;U2 + 2 (Ai + A3) u^) + 2 u (cTi + 2 Re(<T2) + 0-2) + 6 vq Re{as) sin(3 a) 



+ 2(4Re(Kl+K2) + 2(Kl + K2)+T2)tif^ 



+ 2vov^ 



+ 2uv^ 



3in(Q;) (2-3 sin2(/3)) - C03(q;) sin (/3) 003(27) 



yi 



■ cos (2 a) (2-3 sin^(/?)) + Vs sin (2 a) sin2(/?) 003(27) 



y2 



In the following sections we present the Higgs mass matrices Ai^ for the two minima 
around which we have perturbed to find our numerical solutions shown above. We use the 
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following parameterization for the SM-like Higgs doublets 0: 



VEV + 0'' + i 0* 
0^'' + i 0" 



We define M"^ as: 
'''' 



all fields = 



where 0, G {00,01,02,^1,6,6} and x e {r,i,cr,ci}. 

We give the mass matrices in the basis {^f, , , 0f , 0f , 0o} where — r,i,cr,ci for the 
different components of the Higgs doublet fields. For the calculation of the mass matrices 
we have assumed that all the parameters in the Higgs potential are real such that there 
is no mixing between the real and the imaginary parts of the components of the Higgs 
doublet fields. 



B.2 Mass Spectrum for the Minimum (^j) = (01,2) 

The non-trivial VEV conditions in this case are: 

2 

v{2V3vvoRe{K6) + -{3Xl+4Xl)v'^ - nl + {2Re{K5) + K5 + Ti)v^) = 



= and 



V3 



Re(«;6) + 2XoVq - vq + v'^ vq (2Re(«:5) + + ti) 







The resulting mass matrices for the Higgs fields have the following structure: 



( nil 'm2 ni2 
mi m2 \/3m 
mi — -v/Sn 



V 



2 me 


ma \ 


me 


ms 


me 


ms 








m4 







m5 / 



The eigenvalues of such a matrix are given by: 

^ ^mi + 2 m2 + ms ± ^ (mi + 2 m2 — m^Y + 12 m|^ 

- ^mi — m2 + m4 it ^ (mi — m2 — m4)^ + 24 mg^ each two times 
The corresponding characteristic polynomial is given by: 

[-6 ml + (mi - m2 - x) ("i4 - x)] ^ [-3 m| + (mi + 2 m2 - x) ("is - x)] =0 
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For the fields 0'^'' and 0'^* the variables rrii with i — 1, 6 have the following form: 



mi = 2{2v'^ Xf- 1^1 + VqTi) 
1712 = -f (v3foRc(K6) + 2A|f) 

ma = -V {VSvo (2 Re(K5) + K5) + 2 Re(K6) 

O 

7714 = 2 (— + ^0 '^1 + ''"2 v"^) 

m5 = 2 (2Aoi;o - /if + rii;^) 

For the uncharged scalar fields 4>'^: 
2 

mi = 2(-v^(5Af + 8A| + 4A|)-/xi + vg(2Re(K5) + K5 + ri)) 

m2 = -V (3 V3 Re(/«6) uo + 2 (Af - 2 a| + 2 a|)) 
o 

'TIS = ^ V (3Re(K6)'u + \/3 (2Re(«;5) + k;5 + Ti) i;o) 

m4 = 2 {-i4 + (o-i + 2 Re((T2) + c^2) ^^o + (4 Re(K;i + K2) + 2 {ki + k2) + T2) v'^) 

m5 = 2 (6 Ao^o ~ A*i + (2Re(K5) + /^s 

4 

me = ^ fQf Re(u;i - 0^2 — wa) 
v3 

and for the uncharged pseudo-scalars 0*: 

mi = ^v2(3Af -4A|) -2//| + 2i;2(k5 + ti -2Re(At5)) 

m2 = ^v{V3Re{Ke)vo + 2{Xl + Xl)v) 

4 ^ 
= (2V3foRe(K5) + Re(K6)t') 

m4 = 2(-/i2 + (cTi - 2Re(cr2) +a-2)?;o - (4Rc(ki + K2) - 2(k;i + K2) - 7-2)?;^) 
mg = 2(2AoWo-/^i + ('^5 + ri-2Re(K5))t;2) 

me = — - V (Re(K3 - •\/3k4) + V3Re(u;2) '"o) 



B.3 Mass Spectrum for the Minimum (0o) = ^o? (^2) = 

(6> = ^ and (6,3> = (0i> = 
The three non-trivial VEV conditions are: 

V (2 v"^ {X{ + AXl)-nl + vl (2 Re(K5) + K5 + n) + (4 (2 Re(Ki) + ki) + r2 + 2 rg) 

+AuvQ'Re{ijj2+ oj^ — oJi)) = 
2 (Ai + A3) - 3 vo Re((T3) - /xi u + u (ai + 2 Re((T2) + ^2) 

(4(2Re(Ki) + Ki) + r2 + 2r3) + 2i7^i;oRe(a;2 +a;3 - wi) = 
2 Ao - vq + Vo {(Ti + 2 Re(cr2) + ^2) - Re(a3) + 1/ vq (2 Re(At5) + K5 + n) 

+2t>^ 'uRe(a;2 + a;3 — wi) = 
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The mass matrices for the Higgs scalars have a block structure: 



( nil mg ^ 
m2 ms 
. m2 
. m4 

1715 1^7 

\ me J 

The corresponding eigenvalues are: 

1714 , 'Ti2 i m3 

and the solutions of the characteristic polynomial: 



mi — X ms mg 
mg m-j ruQ-x 







For the fields (j)'^'^ and (j)" the variables rrii with i = 1, 9 have the following form: 

mi = 2{2v^ {\\ + A \l) - 1x1 + vln+u^ {t2 + 2 T^i) - 'iuvQKQ{uJi)) 

m2 = 2{2v'^ {\\-2\l)- ^il+vlTl+U^{T2-T^)+2uVoKQ{ljOl)) 

1713 = 4t; (2 -u Rc(k3) + fo Re(«;6)) 

m4 = 2(2«2(Ai-A3)-;u| + vgcJi+v2(r2-2r3) + 2Mt;oRe(o-3)) 

1715 = 2{2u^ + Xs) - nl + v^ai + v^ iT2 + 2T3) -2uvoRe{a3)) 

niQ = 2{2vlXo- III +u'^ai+v'^ Ti) 

m-r = 2{uvo{2Re{a2) + a2)-u'^Re{a3)-2v^Re{LOi)) 

1718 = Av{2u{2Re{Ki) + Ri)+voRe{LLi2+L03)) 
mg = 2v{vo{2Re{K5) + kq) + 2uRe{uj2+L03)) 

For the uncharged scalar fields 0^: 

mi = 2 (6 (Af + 4 A^) - /ii + (2 Re(/i5) + /i5 + n) + (4 (2 Re(Ki) + ki) + r2 + 2 T3) 

+4:UVo Re{uj2 + W3 - uJi)) 
1712 = 2{2v'^ {Xl- 2X1 + 2X1) - fil + v^ {2 Rc{k5) + kr^ + Ti) + u'^ {2 Reim) 

+Ki + 3(2Re(«;2) + ^2) + T2 - T3) + 2uvoRe{LJi - UI2 - 0J3)) 
ma = 12vvQRe{K%) 

rriA = 2 {2u^ (Ai + A3) - /X2 + (^^i + 2Re(a2) + 02) + (4 (2Re(K2) + ^2) + r2 - 2T3) 
+Quvq Re(cr3)) 

m5 = 2 (6u^ (Ai + A3) - fil + vl ((71 + 2Re(o-2) + ^2) + (4 (2Re(Ki) + ki) + T2 + 2T3) 
Re(cr3)) 

me = 2 (6 vl Aq - /x? + (cti + 2 Re((72) + ^2) + (2 Re(K5) + K5 + n)) 
mr = 2(2v2Re(cj2 +^3 - wi) + 2«i;o (o-i + 2Re((72) + 0-2) - 3'u^Re(cr3)) 

= 4t; (n(4(2Re(Ki) + Ki) + r2 + 2r3) + 2?;oRe(a;2 + a;3 - wi)) 
mg = 4z; (vo (2Re(K5) + + Ti) + 2uRe(u;2 + <^3 — <^i)) 
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and for the uncharged pseudo-scalars 0': 

mi = 2 (2 ^2 (Af + 4Xl)-nl+ (ks - 2 Re(K5) + n) + (4 (aci - 2 Re(Ki)) + T2 + 2 rg) 

—4 u vq Re(a;i +L02 — co^)) 
1712 = -2{2v'^{-X{ + 2Xl + 2Xl)+f4 + VQ{2Re{K5)-K5-Ti)+u^{2Re{Ki) 

—ki + 3 (2Re(/«2) — K2) — T2 + T3) — 2ui;oRe(a;i + UJ2 — UJ3)) 
ma = 4v {2uKe{K2 — Vs K4) + VQ'Re{KQ)) 

1714 = 2 (2 (Ai - 2 A2 - A3) - /ii + (f^i - 2 Re((72) + ^2) + v'^ (4 (k2 - 2 Re(K2)) + r2 - 2 rg) 

+2nt;o Rc(cr3)) 

ms = 2(2M2(Ai + A3)-/xi + u^(ai-2Re(c72)+cr2)+i;^(4(Ki-2Re(Ki))+T2 + 2T3) 

-2ni;o Re(c73)) 

me = 2(2t)gAo-^? + ti^(cri-2Rc(o-2)+cr2)+f^(K5-2Rc(K5) + Ti)) 
my = — 2 (u^ Re((j3) + 2^;^ Rc(a;i + a;2 — '^3) — 4mvo Re(c72)) 
mg = S-y (4'uRe(Ki) + fo Re(u;2)) 
mg = (uo Re(K5) + 'uRe(a;2)) 
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